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The evolution of the Cabibbo-Kobayashi-Maskawa matrix and the quark Yukawa couplings is 
performed for the one-loop renormalization group equations in the universal extra dimension model. 
It is found that the evolution of mixing angles and the CP violation measure J may rapidly vary 
in the presence of the Kaluza-Klein modes, and this variation becomes dramatic as the energy 
approaches the unification scale. 



PACS numbers: 12.15.Ff, ll.10.Hi, ll.10.Kk 



I. INTRODUCTION 

With the Large Hadron Collider (LHC) now up and running, exploration of the realm of new physics that may 
operate at the TeV scale has begun. Among these models those with extra spatial dimensions might be revealed in 
such higher energy collider experiments, where the universal extra dimension (UED) model makes for an interesting 
TeV scale physics scenario; as it features a tower of Kaluza-Klein (KK) states for each of the standard model (SM) 
fields, all of which have full access to the extended spacetime manifold This particular scenario has recently been 
extensively studied in the literature, such as investigations of electroweak symmetry breaking, proton stability, gauge 
hierarchy and fermion mass hierarchy problems, B physics, dark matter etc. [2|-|7|. This model has been a fruitful 
playground for addressing a variety of puzzles in the SM. 

On the other hand, it is well known that in the SM, the quark sector's flavor mixing is parameterized by the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix: 

V ud V us V ub \ fV u V V2 V 13 
Vckm = | Vcd Vcs Vob = V 21 V 22 V 23 I , (1) 

v td v u v tb J \ v 31 v 32 v 33 

which makes it possible to explain all flavor changing weak decay process and CP-violating phenomena to date, where 
the 10 year run of Babar at SLAC and the Belle detector at KEK has greatly improved our knowledge of the CKM 
matrix elements. In particular, for the standard parameterization of the CKM matrix, which has the form: 

C12C13 S12C13 s 13 e~ tS 

V C km = \ si 2 c 23 - c 12 s 23 s 13 e lS c 12 c 23 - s 12 s 23 s 13 e tS s 23 c 13 (2) 
S12S23 ~ ci 2 c 23 si 3 e lS -c 12 s 23 - s 12 c 23 si 3 e l& c 23 c 13 

where s\ 2 = sin#i2, c\ 2 = cos6\ 2 etc. are the sines and cosines of the three mixing angles #12, 23 and #13, and 5 is 
the CP violating phase. 

The CKM matrix highlights that, in the quark sector of the SM, we have ten experimentally measurable parameters, 
i.e. six quark masses, three mixing angles, and one phase. A completely satisfactory theory of fermion masses and 
the related problem of mixing angles is certainly lacking at present, however, there has been considerable effort to 
understand the hierarchies of these mixing angles and fermion masses in terms of the renormalization group equations 
(RGE) fcj-tlij. Recall that in order to explore the physics at a high energy scale we use RGE as a probe to study 
the momentum dependence of the Yukawa couplings, the gauge couplings, and the CKM matrix elements themselves. 
As such we can consider one of the primary goals of the LHC as being to uncover any new dynamics within the 
TeV range. Instead of assuming the RGE goes from the Mz scale up to the GUT scale (10 15 GeV) by using the 
SUc{3) x SUl(2) x C/y(l) symmetry, we know that models with extra dimensions may bring down the unification to 
a much lower energy scale. 



'Email: alan.corncll@wits.ac.za 
^Email: luxin.liu9@gmail.com] 



2 



The UED model we shall consider places all SM particles in the bulk of one or more compactified extra dimensions. 
In the simplest case, there is a single flat extra dimension of size R, compactified on an Si/Z 2 orbifold. Therefore, 
from a 4-dimcnsional view point, every field will have an infinite tower of KK modes, with the zero modes being 
identified as the SM states. If these KK modes are indeed within the TeV range, they would modify the running of 
the RGE at relative low energy scales. However, like any higher dimensional theory, the UED model should be treated 
only as an effective theory which is valid up to some scale A, at which a new physics theory emerges. Between the 
scale R^ 1 where the first KK states are excited and the cutoff scale A, there are finite quantum corrections to the 
Yukawa and gauge couplings from the AR number of KK states. Up to the scale the first step KK excitation 
occurs, the RG evolution is logarithmic, controlled by the SM beta functions. With the increasing of the energy, that 
is, when the KK threshold is crossed for each successive mode, new excitations come into play and govern new sets of 
beta functions. The values of physical parameters such as Yukawa couplings and gauge couplings do not run in the 
old SM fashion, instead they receive finite quantum corrections whose magnitudes depend explicitly on the value of 
this cutoff parameter. As a result, once the KK states are excited, these couplings exhibit power law dependences on 
A. This can be illustrated if AR > 1, to a very good accuracy, the generic SM beta function is shown to have the 
power law evolution behavior Q 

pSM ^ pSM + (5(/i) _ x) p f (3) 

where /3 is a generic contribution from a single KK level, and where its coefficient is not a constant but instead S(ii) = 
/j,R, with n Max = A, reflecting the power law running behavior. As a result of faster running, the gauge couplings 
tend to lower the unification scale down to a relatively low order, which might be accessible to collider experiments, 
such as the Tevatron Run-II, LHC or the proposed international linear collider (ILC). Therefore, constraints from 
precision electroweak tests and current (or future) collider data would yield bounds on the compactification radius R. 

The RGE are an important tool for the search of the properties of the quark masses and the CKM matrix at 
different energy scales. It is therefore of great interest to have an implementation of the UED model in studying these 
RGE. In this paper, we consider the UED model with a single compactified extra dimension. In section 2, we start 
from the KK expansion of the SM fields and show the contributions of the relevant one-loop diagrams to the Yukawa 
couplings at each KK level. The RGE of the Yukawa couplings are then derived by using the anomalous dimensions 
of the wave function renormalization of the relevant fields and the vertex renormalization. We shall then derive the 
one-loop RGE for the CKM matrix. In section 3 we shall quantitatively analyze the evolution of these RGE from low 
energies up to the unification scale. The scale dependence of the mixing angles as well as the CP violation measure 
J will be plotted. We also calculate their evolution behaviors for different compactification radii R. The last section 
is devoted to a summary and our conclusions. 



II. RENORMALIZATION GROUP EQUATIONS 

Starting with the SM, which is based on the group structure SUc{3) X SUl(2) x E/y(l) with one Higgs doublet, 
the mass matrices arise from the Yukawa sector of the theory as given by 

C Y = Y v eHuQ + Y D H*dQ + Y E H*eL , (4) 

where Q and L are the SUl{2) doublets for the quark and lepton sectors, u, d and e are the right-handed SUl{2) 
singlets for up-type, down-type quarks and leptons respectively, H is the Higgs doublet, and e is the 2x2 antisymmetric 
tensor, with £12 = —£21 = — 1- The above equation implicitly contains the summation of the quark generations as 
well as the summation over the SUl(2) indices. The information of the physical observables at the scale Mz can be 
extrapolated to a higher energy scale by means of the RGE. It is well known that the evolution of the generic Yukawa 
coupling 

Y$rM , (5) 

which describes the fermion-boson interactions, is given by the beta function. Although the bare constants are 
independent of the renormalization scale, the renormalized coupling constants will depend on the choice of the scale 
parameter /1. As a result, the Yukawa coupling renormalization depends on the corresponding beta functions, including 
contributions from the anomalous dimensions of the field operators. That is, its evolution is given by: 

^ ln yR = 5% ln Zh + l fi IH ln z ^ + l^lk ln z * ~ ln Zcoupin9 ' (6) 

where Y R is the renormalized Yukawa coupling constant (we shall drop the index R for the remainder of the paper), 
and Z^ L , Z^ R and Z§ are the wave function renormalization constants related to left-handed, right-handed fcrmions 
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and Higgs boson respectively, and Z coup u ng is the vertex rcnormalization constant; or in terms of the anomalous 

Id d 
dimensions, Jwave = -^p-^-^Z wave ,j coup ii nci = U-^-lnZ COUP ii n „. In the current context, we have chosen to work with 
2 Ofi op 

the minimal UED model, i.e. the extra dimension is compactified on a circle of radius R with a Z 2 orbifolding, which 
identifies the fifth coordinate y — > — y. The 5-dimcnsional KK expansions of the weak doublet and singlet as well as 
the Higgs field are shown (the corresponding coupling constants among the KK modes are simply equal to the SM 

Y 5 

couplings up to normalization factors, e.g. Yjj — u ) as below: 



'■kR 



H{x,y) 
u(x,y) 

Q{x,y) 
d{x,y) 
L(x,y) 
e(x,y) 



1 



'■kR 
1 



'txR 
1 



H(x) + V2^1/„(x)co S p|)l , 

n=l ) 

oc 

u R (x) + V2Y, [u n n(x) cos (^) + sin 



71=1 



'■kR 
1 



'txR 
1 



q L (x) + [QIW cos (^) + Q n R (x) sin 

k n— 1 

OO 

d fl (x) + y2^ [d£(s)a*(I|) +d2(x)sm(!| 



71 = 1 
OO 



1 



L L {x) + [Ll(x) cos (^) + L R {x) sin 



n=l 



'■kR 



OO 



X cos 



+ e£ (cc) sin 



(VOL 



(7) 



The zero modes in the above equations are identified with the 4-dimensional SM fields, whilst the complex scalar field 
H is Zi even field, and there is a left-handed and a right-handed KK mode for each SM chiral fcrmion. Note that 
in models with UED momentum conservation in the extra dimensions, we are led to the conservation of KK number 
at each vertex in the interactions of the 4-dimensional effective theory (or strictly speaking, the KK parity (—1)" is 
what remains conserved, where n is the KK number). In the bulk we have the fermion and gauge field interactions 
as follows: 

kR 



^Leptons = / dy{iL(x,y)T V M L{x,y) + ie(x,y)T V M e(x,y)} 



o 

jvTI 



C 



Quarks 



dy{iQ(x, y)T M V M Q(x, y) + iu(x, y)T M V M u(x, y) + id{x, y)T M V M d(x, y)} 



Mr 



(8) 



where T = (7^, 17 ), and M = 0, 1, 2, 3, 5. Explicitly, the kinetic terms are given by: 

DMQ(x,y) -- 
V M u{x,y) -- 
V M d(x,y) -- 



d M + ig%G M + ig%W M + i g3i#A/ J Q{x, y) , 

9m + iglG M + i^g\B K ^j u(x, y) , 
9m + iglG M + i-^-gfB M j d(x, y) , 



dM + iglW M +i— g\B M 



V M L(x,y) -- 
D M e(x, y) = (d M - ig\B M ) e(x, y) . 



(9) 



The gauge couplings 33, g\ and g\ refer to those of the SU(3), SU{2) and U(l) gauge groups respectively, and are 



related to the 4-dimensional SM coupling constants gi 



<ttR 



After integrating out the compactified dimension, 



FIG. 1: The one-loop corrections of the additional diagrams (to the SM type diagrams) from the fifth component of the vector 
fields to the Yukawa couplings in Eq.^jj), introduced at each KK excited level. The dashed line is for the Higgs field, the dotted 
line is for the A$ scalar. In TaWeQJfche contributions to the anomalous dimension for each these possible diagrams is presented. 



TABLE I: The anomalous dimensions for each diagram in FigQ] JVote that the columns refer to each type of fifth component 
for the vector fields (pi , gi and gs), for each type of wavefunction renormalization and proper Yukawa vertex renormalization, 
along with their corresponding contributions in the rows for each type of Yukawa beta function. 









'Coupling 


Yukawa beta function 


9i 


92 


53 


9i 


52 


53 


91 


92 


53 


Up — type quark 
Down — type quark 
Lepton sector 
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2.91 5 


1„ 23 
292 3 

1 _ 23 
2.92 3 
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92 5 
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253 g 
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253 g 


1 _ 2/2N2 3 
251 (3) 5 
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251 (-3) 5 
1-23 

s5i 5 




1 „ 2 8 
253 g 

1 „ 28 
253 e 


a 2 2 3 
51 9 5 
. 213 
"51 95 

r, 23 

51 5 




r, 28 

53 3 

r, 28 

53 3 



the 4-dimcnsional effective Lagrangian has interactions involving the zero mode and the KK modes. However, these 
KK modes cannot affect electroweak process at tree level, and only contribute to higher order electroweak processes. 
For the one-loop diagrams of the Yukawa couplings, we choose the Landau gauge in what follows, as many one-loop 
diagrams are finite in the Landau gauge and have no contribution to the renormalization of the Yukawa couplings. 
We therefore consider the RGE for the quark-Higgs Yukawa couplings from which we obtain the evolution of the 
quark masses and the CKM matrix. The one-loop Feynman diagram contributions to the Yukawa couplings in the 
SM and UED model have been explicitly illustrated in Ref.[l, H|. In the UED model, where for each energy level n,, 
we effectively have a heavier duplicate copy of the entire SM particle content. That is, at each KK excited level, the 
KK tower corresponding to the fields in Eq.([7]) exactly mirror the SM field ground states. However, new contributions 
from the A5, 

As(a; ' y) = v^^^ (x)shi (D ' (10) 

n— 1 

interactions (that of the fifth component of the vector fields, i.e. ^5 = {G5, W5, -B5}), as illustrated in FigfTJ also 
contribute. In contrast, the fifth component of the gauge bosons Ar,(x,y) is a real scalar and does not have any zero 
mode, transforming in the adjoint representation of the gauge group. 

That is, each of the other graphs which contribute, due to the new excitations, exactly mirror the zero mode SM 
ground state, where their contributions to the anomalous dimensions are exactly the same as those in the usual 
4-dimensional SM. However, there is a subtlety here, due to the existence of both left- and right-handed chiral KK 
modes in Eq.([7]) (when the expansion of fermion fields is done), that is, we need to double count the anomalous 
dimension contributions in contrast with those of the SM zero modes. We use dimensional regularization to extract 
the anomalous dimensions from the divergent parts of the above one-loop diagrams, where in TableUwe list the results 
of the anomalous dimensions for the relevant fields and vertices in Fig[T] related to the up-type, down-type quarks, 
and lepton Yukawa couplings respectively; where for simplicity we have omitted a common multiplicative factor of 

k in the table. The one-loop RGE for the Yukawa couplings take the following form: 

16. 2 ^ = ft" + /%**>, 

I6n^ = (11) 
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where the Yukawa coupling beta functions have the form: 



P S u M 



oUED 

Pu 



792 



11 



< (S(t) 
+2{S{t) 



17 

" 20 
28 



91 



9 

93 



3 

" 2 

15 2 
Y 9 2 



l)Tr 



101 



120 

ZYlY D 



91 



Y ] d Yd 
3 

+ 2 
YtY E 



Tr 



}■ 



Y„ 



ZY^jYu 



YlY D 



ZYlY D 



YlY E 



(12) 



and 



oUED 

Pd 



= Y E 



89! + \gl 



{ (S(t) - 1) 

+2(S(t) - l)Tr 



1 2 
Jffi 



28 2 

3 y3 



Y, 



15 2 

^52 



17 
120 



3rlY D 



Y E Y E 



Tr 



ZYIjYu 



YlY D - Y^Yjj 



3Y^Yd 



Y e Ye 



(13) 



Note that here t = \n(fj,/Mz) is the energy scale parameter, and S(t) = e l MzR, where we have chosen the Z 
boson mass as the renormalization point. The coupling constant g\ is also chosen to follow the conventional SU(5) 
normalization. 

Note that in deriving Eqs. (|12ll3[) . for the factor of g\ in the UED beta function of Eq. (fT2")) . as an explicit example, 
we have 7^ L + 7^ — ^coupling = —4/3^3 from Table Q] Together with a factor — 8c/| which is read off from the first 
line of Eq. p^| . This results from other one-loop graphs of the KK modes which mirror those of the SM zero mode. 
This then leads to a total factor — 28f/|/3. At this point, our results differ from those of Ref. [1]. In Table fl] we 
have positive signs for the wave function anomalous dimensions. Conventionally, the sign of ^ wa ve can be justified 
by considering the fact that the wave function renormalization constant is less than unity for the gauge independent 
interacting field [l5j], and the divergent parts of these one-loop diagrams are independent of A 5 gauge dependencies. 
Note that while the zero mode fcrmions are chiral as a result of orbifolding, the KK quarks and leptons at given levels 
are vector like. This accounts for the relative factor of 2 between the first and second terms of the UED beta function 
in Eqs. (|12ll3[) for the overall proportionality factor 5—1, since both the KK left- and right-handed chiral states can 
simultaneously contribute to the closed fermion one-loop diagrams. The interaction between A§ and the Higgs field 
only contributes to the renormalization of the Higgs mass, thus the wave function anomalous dimensions of the Higgs 
field is immune to the effects of A§. Similarly, for the lepton sector we have: 



I6n 2 



,dYE 
dt 



P. 



SM , oUED 

' Pe 



(14) 



in which 



oSM 

Pe 



nUED 

'E 



Yt, 



-92 



-91 



Ye (S(t) - 1) 



-YIYe 



Tr 



15 



9% 



99 
40 



91 



ZYlYu 

\y e y e 



Y E Y E 



2(S(t) - l)Tr 



ZY T \Y V + SYXYd 



Y e Ye 



(15) 



Or, explicitly in diagonal form 



16tt 2 ^ = y a 2 [2(25(t) - 1)T - 2G E + 3S(t)y a 2 } 

/9 9 \ /15 99 

where we use Y E = diag(y e , y M , y T ), G E = f -gl + -g 2 j + (S{t) - 1) f —g 2 + —g\ 

Y E Y E ]. Eqs. (|llll6|) constitute a complete set of coupled differential equations for the three families. 



(16) 



and T = Tr[3Y^Yu + 3YlY D 



V, 



Further, after diagonalizing the square of the quark Yukawa coupling matrices by using two unitary matrices U and 



UYfruUi = di ag (f 2 J 2 J 2 ) 
VY^Y D V^ = diag (h 2 ,h 2 ,h 2 ) 



(17) 
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in which f 2 , f 2 , f 2 and h 2 d , h 2 , hi are the eigenvalues of YyYjj and Y^Yjj respectively, it follows that the CKM matrix 
describing the quark flavor mixing in the charged current is given by 



VCKM = UV^ 



(18) 



We will now obtain from Eq. ([TTj) the RGE of the elements of the CKM matrix. Note that we can see from these 
equations that one needs to know the running of Eq. to obtain the evolution of the CKM matrix. By imposing 
the unitary transformation, Eq. (|17p . on both sides of the evolution equations of Y^Yjj and Y^Yn, and taking the 
diagonal elements, we obtain the following two relations: 



16tt 



dt 



= ff 



16ir 2 



> dhj 



dt 



2(2S(t) - 1)T - 2Gu + 3S(t)f 2 - 3S(t) E h 2 \V l3 

3 

2(2S(t) - 1)T - 2G D + 3S(t)h 3 2 - 3S(t) £ f 2 \V l3 



(19) 



which describe the evolution equations of f 2 , f 2 , f 2 and h\, h 2 sl h\ respectively, where 



Gjj 
Gd 



8^ + |ffS + ^S? + (5(t)-l) 



28 



-5 3 



15 



9l 



91 



8ff 3 2 



101 
120 

+ ^3i 



1 



28 



15 



(20) 



Recall that the eigenvalues of the Yukawa coupling matrices are the ratios of the fermion masses to the Higgs vacuum 

v 

exception value v, i.e. m = —j=Y (where v — 246 GcV). The quark mass matrices then appear after the spontaneous 

v2 

symmetry breaking from the quark-Higgs Yukawa couplings, where Eq. (|19p describes how the squares of the physical 
Yukawa couplings evolve. Considering the variation of the square of the Yukawa couplings, we may impose two new 
unitary matrices to make them diagonal. Thus, by applying Eq. (|18[) . we are led to the variation of the CKM matrix 
and thus its evolution equation (beyond the threshold i? _1 ): 



16tH 



>dVu 
dt 



-s(t) 



E fi 2 + fj 2 , 2 T/ v * v , sr^ 
2 r 2 m Vlm V 3 m V 3 k ' 
Ji ~ J j 



j,m=£k 



hk 2 + h m 2 2ta 

,2 i 2 J 3 v 3m v jk Vim 
Ilk ~ "-m 



(21) 



The RGE for the squares of the absolute values of the CKM matrix elements, i.e. the rephasing invariant variables, 
can now be calculated, where one can easily show the following relation holds: 



16^1^ 



dt 



S(t) { 3 \V t3 1 2 iff + h 2 - ]T f% \V k3 \ 2 -J2 hi \V>k\ 2 )-3f 2 J2 



1 



k^i 



ft-fl 



■{2h 2 \V k3 \ 2 \V l3 \ 2 



E^w 



-3^ 2 E 



h 2 



2/ 2 l^l 2 l^l 2 + EA 2 ^ 

k^i 



(22) 



where 



Vi 



iklj = 1 - \V tl \ 2 - \V H \ 2 - \V k3 \ 2 ~ \V l3 \ 2 + \V a \ 2 \V k3 \ 2 + \V kl \ 2 \V i3 \ 2 



(23) 



Note that this evolution does not explicitly depend on the gauge couplings, but on the evolution of the Yukawa 
couplings. 

The structure of the one- loop RGE for the gauge couplings is given by: 1 



16^ = 
dt 



SM 



+ (S(t) - 1% 



9i 3 , 



(24) 



1 Here our calculations agree with the results of Ref. since we need to double count the contributions from both the right- and 
left-handed KK modes for the internal fermion loop, which is in contrast with the results of Ref. 0] and [l6| . 
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where bf M = ^^9' — 6~~'~ an< ^ ^ = (lo'6"' — 2~) *° ^ e coll * r ^ u *' ons n ' om the KK modes, the gauge 

couplings are also altered to a power law evolution. Since the appearance of the extra dimensions causes their running 
to vary much more rapidly, it may well be possible to contemplate a scenario in which perturbative gauge coupling 
unification is preserved to a much lower energy scale. Therefore, if the unification scale is sufficiently low we can 
actually avoid the gauge hierarchy problem. 



III. PROPERTIES OF RGE EVOLUTION 



From this full set of one-loop coupled RGE for the Yukawa couplings and the CKM matrix, together with those 
for the gauge coupling equations, one can obtain the renormalization group flow of all observables related to up- and 
down-quark masses and the CKM matrix elements. We assume the fundamental scale is not far from the range of 
LHC scale, and set the compactification radii to be Rr 1 = 1 TeV, 2 TeV, and 10 TeV respectively. In the limit when 
the energy scale is much smaller than since the energy of the system is less than the excitations of the first KK 
modes, the theory reduces to the usual 4-dimensional SM, and the existence of the KK excitations are ignored. When 
/j > Br 1 , excitations of many KK modes become possible, and the contributions of these KK states must be included 
in all physical calculations. This is characterized by the second term in Eq.Q in the general beta function. Hence we 
apply the full RGE Eqs. (|16ll9l22j) to explore the scaling dependence behaviors of these physical observables. For the 
gauge couplings we take the initial inputs a\(Mz) = 0.01696, a 2 (Mz) = 0.03377, and az{Mz) = 0.1184. Once the 
energy passes the excited KK modes tend to increase rapidly this running of the gauge couplings, and ultimately 
change the scale dependence of the gauge couplings from logarithmic to those of a power law as a function of /i. 
Quantitatively, due to the fast running of the gauge couplings, we find they nearly meet at around t = 5.8, 6.5, 8.2 
(that is, for fi ~ 30, 60, 330 TeV) for radii R~ l = 1, 2, 10 TeV respectively. 

The extra dimensions naturally lead to gauge coupling unification at an intermediate mass scale. To illustrate the 
power law dependence of the Yukawa couplings quantitatively, we take m u (Mz) = 1-27 MeV, m c (Mz) = 0.619 GeV, 
m t (M z ) = 171.7 GeV, m d (M z ) = 2.90 MeV, m s (M z ) = 55 Mev , m b (M z ) = 2.89 GeV, m e (M z ) = 0.48657 MeV, 
Tn^Mz) = 102.718 MeV, and m T {Mz) = 1746.24 MeV [T3| for the initial input values for fermion masses, and run 
the RGE of the Yukawa couplings from Mz up to the GUT scale for our three different compactification radii Rr 1 = 1 
TeV, 2TcV, and 10 TeV . The initial Yukawa couplings are given by the ratios of the fermion masses to the Higgs 
vacuum expectation value. 

As illustrated in FigslHEl the Yukawa couplings evolve in the usual logarithmic fashion when the energy is below 
1 TeV, 2 TeV, and 10 TeV for the three different cases. However, once the first KK threshold is reached, the 
contributions from the KK states become more and more significant. The second terms on the right hand side of 
Eq. tfTTj) depend explicitly on the cutoff A, which have finite one- loop corrections to the beta functions at each massive 
KK excitation level. Therefore, the running of the Yukawa couplings, or more precisely, the one- loop KK corrected 
effective four dimensional Yukawa couplings, begins to deviate from their normal orbits and start to evolve faster and 
faster. For the compactification radius R^ 1 = 1 TeV, the Yukawa couplings evolve faster than the other two, reaching 
its maximum value at the unification scale of around 30 TeV, after that point their evolution will "blow-up" due to 
the faster running of the gauge couplings and new physics would come into play. For the radius R^ 1 = 2 TeV, we 
find similar behavior to the R^ 1 = 1 TeV case, where the "blow-up" scale is not very far from that of 1 TeV case. 
However, for the third choice of radius, since the compactification radius is now much higher than the other two, we 
need more energy to push it further toward its "blow-up" point, which is at a higher unification scale. We also observe 
that the Yukawa couplings arc quickly evolving to zero, however, a satisfactory unification of these seems to still be 
lacking. The first generation f u and hd are driven to the order of 10 -6 , while the / t , the heaviest one, is driven to the 
order of 10 _1 . In Ref . (l8| . however, an additional scalar superfield is introduced which is located at the orbifold fixed 
points and therefore lacks a KK tower. This scalar superfield's contribution to the Higgs wavefunction renormalization 
causes the unification of the Yukawa couplings to become feasible. In the UED scenario, the unification of the Yukawa 
couplings is very desirable due to the fast power law running. This feature thus has the potential to address the 
problem of fermion mass hierarchy. We plan to come back to this point in a future work. As such, we have so far 
observed the Yukawa couplings all decrease with increasing energy, which agrees with what is observed in the SM, 
however, the Yukawa couplings arc driven dramatically towards extremely weak values at a much faster rate. This is 
an interesting feature that distinguishes the UED model from that of the SM. 

We next turn our attention to the quark flavor mixings. Because of the arbitrariness in choice of phases of the 
quark fields, the phases of individual matrix elements of Vckm are not themselves directly observable. We therefore 
use the absolute values of the matrix element \Vij\ as the independent set of rephasing invariant variables. Of the 
nine elements of the CKM matrix, only four of them are independent, which is consistent with the four independent 
variables of the standard parameterization of the CKM matrix. Among these, the complex phase of the CKM matrix 
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FIG. 2: The evolution of the Yukawa coupling fi (i — u top left panel, c top right panel, t bottom panel), where the solid line 
is the SM, the dotted line is the R -1 = 1 TeV UED case, the dotted-dashed line is the 2 TeV UED case and the dashed line is 
the 10 TeV UED case. 



characterizes CP-violating phenomena, which have been unambiguously verified in a number of K — K and B — B 
systems. Conventionally, one uses the Jarlskog rephasing invariant parameter J = lmV u dV cs V* s V* d to present CP 
violation phenomena. Its square can be written as follows: 

J 2 = \Vn\ 2 \V 22 \ 2 \V 12 \ 2 \V 21 \ 2 \Vn\ 2 - \V 21 \ 2 - \V 22 \ 2 - \V 12 \ 2 + \V U \ 2 \V 22 \ 2 + \V 21 \ 2 \V 12 \ 2 ) 2 . (25) 

For definiteness, we choose the \V u b\, \V c b\, \V US \ and J as the four independent parameters of Vckm, and take the 
initial values \V ub \ = 0.00347, \V cb \ = 0.0410, \V US \ = 0.2253 and J = 2.91 x 10" 5 [l9|. In Figg]we plot the energy 
dependence of these four variables from the weak scale all the way up to the unification scale for different values of 
compactification radii R. 

We observe from these plots the following; the CKM matrix elements V u b — 9i 3 e~ lS , V c b — 23 , can be used to 
observe the mixing angles #13 and #23 and that they increase with the energy scale; the variation rate becoming faster 
once the KK threshold is passed. For the mixings related to the third family, the UED effects become sizable and the 
mixing angles #13 and #23 change at a level of 15% between Mz and the unification scale, in contrast with the SM, 
in which the angles only rise by around 5% at similar energy scales. By contrast, the variation of the Cabibbo angle 
appears to be the least sensitive. Due to the smallness of the Yukawa coupling terms on the right hand side of Eg. (|2"2"1) . 
the renormalization group flow of the mixing between the first two families turns out to be very small. Although the 
mixing angle increases all the time, it is rather inert, even in the UED model. It has a maximum variation around the 
order of A 6 , for the Wolfcnstcin parameter A = 0.22. However, for the parameter J, the characteristic parameter for 
the CP non-conservation effects, its variation becomes very significant. The larger the value of the compactification 
radius R, the faster J evolves to reach its maximum. We observe an approximate 30% increase for J at the unification 
scale compared with its initial value. The absolute values of all the remaining magnitudes of the CKM matrix elements 
can be obtained from the unitarity equations, thus we could determine the renormalization group evolutions of the 
full CKM matrix. As depicted in Fig|4] with increasing energy the running of the CKM matrix shows a pronounced 
pattern at the point where the KK modes are excited. As can be seen from Eq. (j2"2"|) . the evolution of the CKM matrix 
is governed by the Yukawa couplings and the factor S(t). They evolve faster in the region where the power law scaling 
of the Yukawa couplings becomes substantial. This effect is explicit for mixings involving the third family, due to the 
large value of their Yukawa couplings. 
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FIG. 3: The evolution of the Yukawa coupling hj (j = d top left panel, s top right panel, b bottom panel), where the solid line 
is the SM, the dotted line is the R~ 1 = 1 TeV UED case, the dotted-dashed line is the 2 TeV UED case and the dashed line is 
the 10 TeV UED case. 



IV. SUMMARY 



UED models with a compactification radius near the TeV scale promise to provide an exciting set of observable 
phenomenology for collider physics. If the compactification radius R is sufficiently large, all SM particles KK partners 
might be detected at accelerators, such as the LHC, which can run at center of mass energies of 10 TeV. In the UED 
scenario we have examined the cumulative contribution of these KK states to the renormalization group evolution of 
the Yukawa couplings and the CKM matrix. Due to the power law running of the gauge couplings, we can bring the 
unification scale down to an explorable range at the LHC scale. This fact is in clear contrast to the SM. We have 
plotted the evolution of the Yukawa couplings in the UED for different compactification radii R, where as depicted in 
Figs[2][3l the rapid decrease of the Yukawa couplings with energy is in clear contrast to the slow logarithmic running 
predicted by the SM. By permitting the Yukawa couplings to receive power law corrections in the UED model we show 
that extra dimensions provide a potential mechanism for explaining the fcrmion mass hierarchy. To introduce new 
fields in the bulk or at the orbifold fixed points that allows new terms in the beta function may give us an alternative 
evolution scenario. A more detailed analysis in this direction will be attempted in a future publication. 

To conclude, we have investigated the consequence of the UED model on the Yukawa couplings and CKM matrix 
elements evolution. We performed a qualitative study of the behavior of different mixing angles and the CP violation 
measure J as well. The energy dependence of | V us | is very weak and qualitatively different from those of the mixing 
behaviors involving the third generation. However, when the unification scale is approached, its variation is rather 
rapid, in clear contrast to the slow variation of the SM. While the evolution of the Cabibbo angle is tiny, the elements 
\V u b\ and \V c b\ increase sizably, the relative deviations for these can be up to 15% in the whole range from Mz to 
the GUT scale. As for the energy scaling of J, the contribution of KK modes is substantial. Its numerical analysis 
shows us that its variation can be raised to more than 30%. The scale deviation of renormalization curves from the 
usual SM one depends closely on the value of the compactified radius R. The smaller the radius is, the higher the 
energy scale we need to differentiate the UED curve from the SM one. A comparison between theoretical predictions 
and experimental measurements will be available once the LHC is running at its full scale and the acquisition of 
data becomes available. It is believed that it will set a strong limit on the parameters of this model, and a precise 
determination of J, \V u b\ or \V c b\ may lead to a discrimination between the SM and extra dimensional models. 
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FIG. 4: The evolution of the CKM matrix element \V u b\i t°P Jeft panel, \V c b\ top right panel, \V ua \ bottom left panel, and J, 
bottom right panel; where the solid line is the SM, the dotted line is the R^ 1 = 1 TeV UED case, the dotted-dashed line is the 
2 TeV UED case and the dashed line is the 10 TeV UED case. 
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